Sensitivity of Nuclear Transition Frequencies to Temporal Variation of the Fine 
Structure Constant or the Strong Interaction 
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There exist in nature a few nuclear isomers with very low (eV) excitation energies, and the combi- 
nation of low energy and narrow width makes them possible candidates for laser-based investigations. 
The best candidate is the lowest-energy excited state known in nuclear physics, the 7.6(5) eV isomer 
of ^^^Th. A recent study suggests that a measurement of the temporal variation of the excitation 
energy of this isomer would have 5-6 orders of magnitude enhanced sensitivity to a variation of 
the fine structure constant {a = 1/137.036) or of a strong interaction parameter {mq/ Aqcd)- We 
reexamine the physics involved in these arguments. By invoking the Feynman-Hellmann Theorem 
we argue that there is no expectation of significantly enhanced sensitivity to a variation in the fine 
structure constant (beyond that obtained from experimental considerations such as the low energy 
and narrow width of the isomer). A similar argument applies to the strong interaction, but eval- 
uating the shift due to temporal variations of the underlying parameters of the strong interaction 
may be beyond current nuclear structure techniques. 

PACS numbers: 23.20.-g,06.20. Jr,27.90.+bb,42.62.Fi 
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The excitation energy of the first excited state of ^^^Th 
has recently been determined [ij to be 7.6±0.5 eV. This is 
the lowest- lying excited state known in nuclear physics. 
Its lifetime is not known, but has been estimated to be 
roughly 5 hours These properties suggest that the 
transition between this state and the ground state of 
^^^Th is within the range of laser sources (in the vac- 
uum ultraviolet) and that properties of the state could be 
probed via high-resolution laser spectroscopy. This pos- 
sibility has led to the intriguing ideaQ] that the transition 
could be used to place stringent constraints on a possible 
temporal variation of the fine structure constant (d) or 
on a temporal variation of a particular strong interaction 
parameter [mq/ Aqcd)- Flambaum^] estimates that a 
temporal variation in the energy of this transition would 
provide 5-6 orders of magnitude enhanced sensitivity to 
a temporal variation in these parameters. The purpose 
of this letter is to examine the physics issues involved in 
temporal changes of nuclear transition frequencies and 
to quantify the corresponding sensitivity to variations of 
the underlying fundamental interactions. 

The energy of a transition (lu) between the ground 
state (Iff.s.)) and an excited state (|/)) of a nucleus is 
determined by the difference of matrix elements of the 
nuclear Hamiltonian, H: 
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where we have introduced the notation ((O)) to denote 
the difference of the matrix elements of the operator O 
between the two states. Any variation in time of the 
transition frequency will depend on the time variation 
of various dynamical constants (denoted generically by 
Xi{t)) embedded in the Hamiltonian, H. Assuming a 
sufficiently slow time variation in the constants, Ai, that 
stationary states are well defined, using the chain rule on 



Eqn. (1) leads immediately to 
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The parameter derivatives of the transition fre- 
quency can be obtained using the Feynman-Hellmann 
Theorem[^ 0], which is sometimes (and appropriately) 
called the parameter theorem. The usual statement of 
the theorem is that if a Hamiltonian H{X) depends on 
a parameter A, the energy shift £'(A) for any (bound) 
eigenstate of that Hamiltonian satisfies Q 



'^'^'-WA)lPi)|»(A)). 



dX 



dX 



(3) 



Since this relation holds for the energy of every (bound) 
eigenstate, it must hold for the difference of two such 
energies. This immediately leads to a very simple, useful, 
and essentially exact form of Eqn. (2) that applies to any 
system 
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For the present problem we split the nuclear Hamilto- 
nian into Hn plus Vc as a convenient shorthand notation 
to represent respectively the strong interaction (plus ki- 
netic energy) and the Coulomb interaction between nu- 
cleons in the nucleus. Let us first restrict ourselves to a 
possible variation only of the fine structure constant. Us- 
ing the fact that the Coulomb potential is linear in that 
constant, we find the very simple result 

^ = ((^c»- = ((/ \Vc\f)- {g.s. I Vc I g.s.)) - , (6) 
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where the variation due to the fine structure constant is 
driven solely by the Coulomb energy difference Q of the 
two states. Forming a fractional frequency shift we 
have 
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where the quantity in brackets can be interpreted as an 
enhancement factor (if larger than 1). 

The conventional description of the two states in ^^^Th 
uses the Nilsson modelfl^. In this picture the states 
are described as an "active" neutron in a deformed well. 
In the present work we diagonalized the 3-D Nilsson 
Hamiltonian for the active neutron using the deforma- 
tion parameters of MoUer and Nix to obtain wave 
functions for the ground state and the isomer of ^^^Th. 
In this model the Coulomb energies of the two states are 
predicted to be the same because the isomer and ground 
state are simply different Nilsson states arising from the 
same deformed nuclear shape. In particular, the pro- 
ton structure of the two states is identical. The Nilsson 
model per se does not incorporate a Coulomb interaction, 
although it is implicitly part of the phenomenological in- 
teraction constants for the active nucleonfiol. [Tl|. Thus, 
in this model ((Vc)) = and there is no sensitivity to a. 

The Nilsson Hamiltonian does not predict these two 
levels of •^■^^Th to be 7.6 eV apart and it is necessary to 
include the first order correction to the Nilsson Hamil- 
tonian (namely the pairing interaction) in order to ex- 
plain the observed splitting of the states in ^^^Th. When 
the pairing interaction is added it plays a crucial role in 
reducing the excitation energy of the isomer of ^^^Th. 
However, {(Vc)) remains zero because t he p roton struc- 
ture of the two states remains the sameflll [T2j. An im- 
proved estimate for ((Vc)) could be obtained by allowing 
the proton deformation of the two states to be different. 
Such a difference would arise if different deformation pa- 
rameters are needed to minimize the ground state and 
the isomer energies, (the Moller and Nix deformation pa- 
rameters were chosen to minimize the ground state en- 
ergy only). Alternatively, new proton-proton or proton- 
neutron terms could be added to the Hamiltonian that 
allow differences in the proton core for the two states. 
The first of these extensions of the Nilsson plus pair- 
ing model is relatively straightforward to implement . 
while the second would require a considerably more so- 
phisticated treatment of deformed heavy nuclei than that 
normally included in Nilsson-like models. The Coulomb 
energy difference will likely not vanish for the real nu- 
cleus and, if experimentally possible, the issue will best 
be resolved by determining the electric quadrupole mo- 
ments of the two states. However, there is presently no 
reason to suggest that ((Vc)) is dramatically larger than 
cj, and in any case the constraint {{Hn + Vc)) = lo must 
be satisfied. 

Extracting the dependence of lo on fundamental con- 
stants in the strong nuclear Hamiltonian Hn is consid- 
erably more complicated. It requires a sufficiently de- 



tailed understanding of that force that the appropriate 
parameters can be identified, and the dependence of the 
Hamiltonian on those parameters must be known in or- 
der to use Eqn. (3). A parameterization in terms of the 
quark mass niq and A^qcd is particularly non-trivial. An 
alternate approach is to consider a parameterization of 
the Hamiltonian in terms of the mass of the nucleon, the 
pion, etc. The simplest parameter relationship is then 
the dependence of the Hamiltonian on the nucleon mass. 
Mm- If we ignore corrections to the nuclear potential 
that are inversely proportional to powers of the nucleon 
mass[l3|, all of that dependence is situated in the ki- 
netic energy 1J|, T . Performing the necessary parameter 
derivatives, one finds 
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Here the ellipsis denotes the dependence on all other 
parameters in the nuclear Hamiltonian. Unfortunately, 
there presently does not exist a model for the structure of 
heavy nuclei that would allow an accurate determination 
of any of these terms, with the possible exception of the 
kinetic energy term. Within the context of the Nilsson 
model the kinetic energy varies from state to state, and 
the scale of this variation is a few MeV. 

Flambaum assumed a simple model for the struc- 
ture of the 5/2+ ground state and the 3/2+ isomer of 
^■^^Th and described the states in terms of the asymptotic 
Nilsson orbits JP[Nn^K] = 5/2+[633] and 3/2+[631]. In 
this model ((T)) = because diagonal matrix elements 
of the kinetic operator only depend on the quantum 
numbers N and and not on A. Hence there would 
be no sensitivity to variations in the nucleon mass in 
this model. However, as discussed above, the structure 
of ^^^Th is considerably more complicated. Using our 
wave functions from the diagonalization of the Nilsson 
Hamiltonian with the deformation parameters obtained 
by Moller and NixfTsJ, we find[16] that ((T)) = -0.49 
MeV. In this model the asymptotic states [633] and [631] 
make up about 45% of the wave functions for the ground 
state and the isomer, respectively. The other important 
orbits contributing to the wave functions are [642], [613], 
[622], [602] and [853] for the ground state, and [642], 
[611], [651], [622], and [871] for the isomer. We find that 
small changes in the deformation parameters from those 
obtained by Moller and Nix lead to significant changes in 
the predicted structure of the two states [l^ because at 
these deformation values several "active" neutron Nils- 
son levels lie very close to one another in energy. These 
changes lead to a correspondingly large change in ((T)), 
and ((r)) is predicted to lie anywhere from zero to a 
few MeV. Thus, even the first term in Eqn. (8) is very 
difficult to predict accurately. 

Reference Q also investigated variation of the strong 
force due to variations of the quark masses, which 
chiral perturbation theory demonstrates is the chiral- 
symmetry-breaking part of the strong force. Although 
the strong force is not proportional to niq in leading- 
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order, pion-mass terms (which depend on niq) are incor- 
porated into the one-pion-exchange potential, which is a 
part of the leading-order nuclear force contribution. Chi- 
ral perturbation theory[l7j provides enough information 
about the nuclear force that estimates could be made of 
appropriate parameter derivatives in Eqn. (3). However, 
an accurate evaluation of Eqn. (4) for such terms is very 
likely beyond the accuracy of present nuclear structure 
techniques. Nevertheless, there are no obvious reasons to 
expect any enhanced sensitivity to temporal variations 
of the strong interaction parameters in to, particularly 
within the framework of the Nilsson model. 

Finally we comment on the origin of the large enhance- 
ment factor in Ref. Similar expressions to ours were 
obtained for the relationship between lu and a (the strong 
interaction case is much more problematic), except that 
a parameter from the Nilsson model (whose value was 
approximately 1 MeV) occurs instead of ((Vc)). There is 
no justification for this. Indeed, if we use the asymptotic 
Nilsson model invoked in or the full Nilsson Hamil- 
tonian plus the pairing interaction, all sensitivity to a 
drops out. 

In conclusion the existence of a huge enhancement fac- 
tor in the fractional frequency shift between the isomer 
and ground state of ^^^Th due to a time-dependent fine- 
structure constant is very difficult to realize in any rea- 
sonable theory of the nucleus. We have examined two 
models for the structure of the ground state and of the 
isomer of ^■^^Th, namely, the asymptotic Nilsson model 
and the full Nilsson Hamiltonian using the deformation 
parameters of MoUer and Nix. In both cases there is no 



sensitivity to a and an extension of the Nilsson model to 
include pairing leads to the same result. The main issue 
is that the magnitude of the enhancement is strongly con- 
strained by the Feynman-Hellmann Theorem. It would 
require two states of nearly identical proton structure 
to have Coulomb energies that differ by 1 MeV. An 
analogous difficulty arises for any enhanced sensitivity 
to variations of the strong interaction parameters. The 
two variations (asymptotic and full) of the Nilsson model 
predict ((T)) to be zero and -0.49 MeV, respectively. Ac- 
curately estimating a frequency shift due to variations of 
strong-interaction parameters, moreover, may be beyond 
current nuclear structure techniques. Of course, the fact 
that u! is very small and that the lifetime of the isomer 
implies a very narrow state would still be important con- 
siderations in selecting the 7.6 eV isomer of ^^^Th for 
laser-based studies. 

We note in closing that many of the arguments pre- 
sented here do not rely on details of nuclear structure. 
They are based solely on the Feynman-Hellmann The- 
orem and thus are equally applicable to studies of the 
sensitivity of transition frequencies in atoms or other sys- 
tems to temporal variations of the underlying parameters 
of the Hamiltonian (such as a or the electron mass). 
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